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摘要 : 方 波 信号 和 双 值 噪声 在 科学 领域 广泛 应 用 ,而 乘 性 信号 在 工程 中 广泛 存在 。 但 是 ， 
就 乘 性 方 波 信号 和 乘 性 双 值 噪声 对 双 稳 系统 的 非 线性 影响 方面 的 研究 还 较 少 ,因此 本 文 研究 
乘 性 方 波 信号 和 乘 性 双 值 噪声 作用 下 双 稳 系统 中 的 随机 共振 现象 。 基 于 绝热 近似 条 件 ,推导 
系统 的 信 噪 比 的 表达 式 并 基于 该 公式 进行 了 不 同 输入 条 件 的 仿真 研究 。 研 究 结果 表明 IER 
比 随 着 乘 性 方 波 信号 和 乘 性 双 值 噪声 幅度 的 变化 可 以 取得 一 个 最 大 值 和 一 个 最 小 值 . 大 的 乘 
性 方 波 信号 幅 值 会 抑制 信 噪 比 , 对 于 噪声 相关 率 4=0.65 情形 , 方 波幅 值 4=1.0 时 的 系统 输 
出 幅度 比 4=0.3 时 的 输出 幅 值 下 降 约 29% ; 而 较 大 的 乘 性 双 值 噪声 强度 可 以 提升 系统 输出 
性 能 ， 当 系统 参数 姑 3.8 , 乘 性 双 值 噪声 强度 为 e=3.5 时 的 系统 输出 幅度 为 E=0.25S 时 的 
输出 幅 值 的 2.5 倍 。 随 着 乘 性 白 噪声 、 加 性 白 噪声 及 系统 参数 的 变化 ， 信 噪 比 曲线 出 现 一 个 
峰值 。 
Ki: 随机 共振 ; 双 稳 系统 ; 乘 性 方 波 信号 ; 乘 性 双 值 噪声 
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随机 共振 是 乘 性 在 随机 动力 系统 中 的 非 线性 现象 。 当 作用 于 系统 的 噪声 、 系 统 输入 信和 号 
和 系统 三 者 间 协 同 作用 时 ， 噪 声 对 系统 输出 性 能 有 所 影响 。 随 机 共振 的 概念 由 文 [1j 提 出 ， 
于 对 周期 性 出 现 的 地 球 冰 川 时 代 进 行 合理 的 建 模 分 析 。 从 此 以 后 , 随机 共振 现象 得 到 广泛 
研究 。 随 机 共振 也 广泛 存在 于 天 文 领域 ， 比 如 ， 太 阳 兆 评 和 喷 流 中 ， 与 磁 重 联 相关 的 随机 
共振 加 速 ， 不 但 产生 电子 ， 而 且 导 致 3He/4He 千 倍 的 增强 外。 小 行星 Bennu (101955) 在 行星 
碰撞 和 随机 共振 的 联合 作用 下 ， 经 过 几 百 万 至 儿 千 万 年 ， 到 达 目 前 的 运行 轨道 ”。 电 子 与 加 
速 圆 环 面 小 区 域 中 的 油 动 磁场 相互 作用 的 随机 共振 加 速 , 可 能 是 人 马 座 A 的 毫米 和 短波 频谱 
的 产生 机 制 ”。 在 射电 天 文 接 收 机 中 也 满足 随机 共振 条 件 ， 即 ， 接收 机 噪声 、 输 入 信和 号 〈 射 
电 天 文 信号 、 无 线 电 干扰 信号 ) 以 及 接收 机 相关 参数 ， 如 增益 、 稳 定性 、 动 态 范围 、 噪 声 系 
数 、 矩 形 系数 等 相互 作用 ， 共 同 作用 决定 了 整个 观测 效果 。 
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方 波 信号 可 以 看 作 在 两 个 电 平 间 切 换 的 数字 信号 。 随 着 集成 电路 的 快速 发 展 ， 方 波 信号 
在 数字 通信 系统 中 起 着 越 来 越 重要 的 作用 。 比 如 ,可 以 作为 数字 多 辑 系统 的 输入 信号 ， 作 为 
时 钟 信号 准确 地 触发 同步 电路 。 另 一 方面 ， 双 值 噪声 是 工程 应 用 中 常见 的 干扰 。 比 如 ， 随 机 
电报 噪声 是 电子 工程 领域 一 种 由 型 的 双 值 噪声 。 在 半导体 器 件 中 ,闪存 产生 的 随机 电报 噪声 
影响 信息 的 正确 读 取 。 

系统 的 输入 信号 通常 具有 加 性 的 形式 ， 即 以 加 法 的 形式 作用 于 系统 。 然 而 ， 在 某 些 情况 
下 ,输入 信号 是 状态 依赖 性 的 ， 输 入 信号 应 该 是 与 系统 状态 变量 相 乘 的 形式 。 目 前 乘 性 信 品 
作用 对 射电 天 文 接收 机 可 能 存在 干扰 特别 是 在 强 干 扰 信 号 作用 下 ， 均 以 乘 性 结果 输出 器 。 
目前 ， 射 电 天 文 台 址 的 无 线 电 环境 测量 较 多 关注 干 护 信 号 频 点 和 品 底 水 平 " ， 但 对 强 干 扩 
可 能 导致 的 接收 机 内 部 交 调 后 的 噪声 变化 却 没有 相关 研究 报道 。 同时 该 研究 结果 ,可 以 用 于 
基于 锁 相 放大 的 微弱 信号 提取 ”中 ， 而 锁 相 放大 后 引起 的 信 品 比 提升 一 直 是 值得 研究 的 问题 
基于 此 , 本 文 骨 在 研究 乘 性 方 波 信号 和 乘 性 双 值 噪声 作用 下 非 线性 双 稳 系统 中 的 非 线性 共振 
行为 。 通 过 仿真 定量 分 析 ， 得 到 非 线性 共振 行为 的 能 谱 分 布 规律 为 强 无 线 电 干 扰 情 况 下 的 
接收 机 设计 、 干 护 信 号 处 理 提供 依据 。 

1 乘 性 方 波 和 乘 性 双 值 噪声 作用 下 的 双 稳 系统 
性 信号 和 乘 性 噪声 作用 下 的 双 稳 系统 满足 au， 
EO la+ FORD) -be + ED) + 7), a) 
sth, a Alb 为 系统 参数 (a>0,b>0)， 
f(t)=r+ec(t)+ As(t) (2) 
HRSS. HO 为 对 称 双 值 噪声 ，G(z)={1,-1} ， 其 相关 率 为 和 s AAWA T HA 


XA 


1, 0<t<T/2 
s(t)= (3) 


=L7T/2<1t<T ` 
r 为 直流 偏 置 6 FA 分别 为 乘 性 双 值 噪声 和 乘 性 方 波 信号 的 幅度 ,车 f(1) =0， 则 系统 (1) 
为 一 双 稳 系统 ， 一 个 非 稳 态 和 两 个 稳 态 x,= 土 Ya/b 。 限 制 参数 f(t) 的 取 值 使 其 满足 
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Fæ a+f(t)>0. EOM 为 相关 白 噪 声 ， 其 均值 为 0， 相关 函数 为 
ii (EDEN, )) = 2DAN -t,) (4) 
{70)7(6)) = 2PEG, =t) (5) 
(E(t,)n(t,)) = (NEG, ))=24V DPS -6) ， (6) 


其 中 ， 刀 和 尸 分 别 为 乘 性 和 加 性 噪声 的 强度 ， 为 两 噪声 间 的 耦合 强度 ，0<4<1 
由 (1)、(5) 和 (6) 式 ， 系 统 (1) 对 应 的 Fokker-Plank 方程 可 以 表示 为 


POD eC, Npa DÀ], me 
ôt Ox Ox 
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F(x,t) =ax—bx* + Dx+ ANDP +xf (t), (8) 


G(x,t) = Dx’ +2AVDPx+P. (9) 


设 方 波 信 号 为 一 个 缓 变 信和 号, 系统 (1) 在 一 个 周期 的 时 间 内 足够 达到 局 域 平衡 , 即 系统 (1) 
满足 绝热 近似 条 件 “， 这 样 可 以 由 (7)- (9) 式 推导 其 长 时 间 的 渐 近 分 布 函数 , 即 


pepe 

a” ACD 
其 中 ，M 为 归 一 化 常数 ， (xn) 为 修正 的 势能 函数 

D60)=| -dx F/G, À). GD) 


exp[—®(x,1)], (10) 


在 绝热 近似 条 件 下 ， 粒 子 在 两 稳 态 的 转移 率 可 以 表示 为 
a 
Wi (1H) = exp[®(x,) — B(x, )] 
+U Jan p 0 + 
a a 
= ex J 22abDP 1 D | . te) 
J2n aE 
a+D+bP/D-—4bA°P/D | P | 
xexp In 
2D Dalb+2AVaDP/b+P 


xexp! # —AVDP -2bAPVDP | D? 
DI —A7)P 
arctan ANP arctan = alb + ANP IND 
DJ(1—A’)P JDU—A’)P 


系统 输出 的 相关 函数 可 以 表示 为 二 
G(t) =[B,(7, £, A) + B,(r, £, A)lexp(—A, t) +B, r, £, De +C(r,€, AEA, (13) 
其 中 


B (r,e, A= [rtetA -mre- A+ Wr+e- A -wr -er a, (14) 
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B,(r,8,A) = Wr +E +A) wr—-e—-A)—-wrt+e—-A)+w(r-e+ A], (15) 
Br, EA) = ZIW +e + Abt re A)-w(r 2 A)— wre ADP, (16) 
CHB, A =AC (HE +A) G(r e-A)+C (r+ 8- AMC (7-64 AY], (17) 


a W (u)-W,(u) 


WO 
SW, CW (4) 
C(u) = * | (19) 
oF WD + WF 
w= ID exp[-i(2j+1)Qt]. (20) 
j=0 


对 (10) 式 两 边 进 行 傅 里 叶 变 换 ， 得 
S(@) =S,(0)+S,(@), (21) 


其 中 
HO LA A)+B,(r,€,A)]+C(r, €, A), (22) 
S,(@)=B,(r, £, A)Q(@) , (23) 
G0)= i+) Ste —(2i +DQ]. (24) 
T i=0 


这 里 ,9(0) 为 零 频 率 处 的 噪声 背景 功率 ; 5,(@) 为 输出 信号 的 功率 谱 。 系 统 输出 信 噪 比 定 义 
为 信号 的 基 波 频率 处 功率 与 背景 噪声 功率 之 比 ， 即 

SNR = 8 B,(r,€,A) | 

m C(r,é,A)+2[B,(r,€, A) + B,(r, €, A)]/4 


(25) 


2 讨论 


基于 绝热 近似 条 件 推导 了 研究 的 双 稳 系统 的 渐 近 分 布 函数 和 稳 态 间 转 移 率 。 由 相关 函数 
= 得 到 了 输出 信号 的 功率 谱 ， 进 而 求 得 了 系统 的 输出 信 品 比 。 现 在 基于 (25) 式 , 分 析 信 品 比 的 
> 非 单调 行为 。 利 用 图 1 和 表 1 AT RE OER FT AA SR LA a). FS 1 和 表 1 可 见 ， 
在 信 噪 比 与 乘 性 噪声 幅度 E 的 关系 曲线 上 出 现 了 两 个 极 值 ， 一 个 极 大 值 ， 一 个 极 小 值 。 随 
ABE 的 增 大 ， 信 品 比 先 快速 增 大 至 最 大 值 ， 然 后 缓慢 地 减 小 到 最 小 值 , 最 后 单调 增 大 。 较 大 
的 乘 性 双 值 噪声 幅度 可 以 提升 系统 输出 性 能 ， 由 表 1 可 得 ， 对 于 3. 8 情形 ，&=3.5 时 系 
统 和 输出 幅度 为 e=0.25 时 的 输出 幅 值 的 2.5 4%; 而 当 4.0 时 ，E=3.5 时 输出 幅度 为 
E=0.25 时 的 输出 幅 值 的 2. 15 倍 。 本 文 研究 发 现 ， 乘 性 双 值 噪声 对 信 噪 比 的 影响 与 文 
[12-14] 中 加 性 双 值 噪声 对 信 噪 比 的 影响 不 同 。 对 于 加 性 噪声 情形 ， 在 信 噪 比 与 双 值 噪声 幅 
度 的 曲线 上 或 者 仅 有 一 个 极 值 ， 即 最 大 值 ”， 或 者 虽然 有 两 个 极 值 出 现 ， 但 先 出 现 最 小 值 
而 后 出 现 最 大 值 ”*。 


-> 


>- 


> 表 1 系统 参数 b 和 双 值 噪声 幅度 = 取 不 同 值 时 的 信 噪 比 (10 ) 
f> Table 1 The SNR for different values of the system parameter b and amplitude ¢ of the multiplicative dichotomous noise 


= E 0. 25 0. 50 1. 00 1.50 2. 00 2. 50 3. 00 3. 50 


í b=3. 8 | 3. 41 5. 22 5. 56 5. 37 5. 19 5. 17 5. 72 8. 52 
b=4. 0 | 4. 39 6. 71 7. 15 6. 89 6. 64 6. 53 6. 95 9. 44 
b=4. 2 | 5. 48 8. 35 8. 93 8. 59 8. 26 8. 06 8. 32 10. 44 
b=4. 5 | 7. 30 11.17 | 11.88 | 11.42 10.94 | 10.58 10. 60 12. 14 
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图 1 系统 参数 b 取 不 同 值 时 信 品 比 与 双 值 噪声 幅度 e 的 关系 曲线 ， 参 数 为 a=5 , D=0.1,P=0.1,4=0.05 ,r=0.2,A=0.1， 
A, =1 


Fig.1 The SNR versus amplitude ¢ of the multiplicative dichotomous noise for a=5 , D=0.1, P=0.1,4=0.05 ,r=0.2,A=0.1, 


A, =1 for different values of the system parameter b 
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由 图 2 和 表 2 研究 信 噪 比 与 乘 性 方 波 信和 号 幅 值 的 依赖 关系 。 由 图 2 和 表 2 可 见 ， 信 噪 比 
随 A 的 增 大 而 非 单 调 变 化 。 随 着 A 的 增 大 ， 信 噪 比 先 减 小 至 一 个 极 小 值 ， 而 后 增 大 至 一 个 
极 大 值 ， 最 后 单调 减 小 。 可 见 ， 相 对 居中 的 幅 值 A (如 对 于 4=0.6S 取 4= 0.65) 可 以 
提升 系统 输出 信和 号, 而 较 高 或 较 低 的 幅 值 A (对 于 多 =0.65 , 取 A>0 或 4=0.3) 则 
会 抑制 系统 性 能 。 比 如 ， 由 表 2 可 得 , 对 于 多 =0.65 情形 ， 信 1. 0 时 系统 输出 幅度 比 全 0. 3 
时 的 输出 幅 值 下 降 约 29%; 24 A=O0.67 IN, 4=1.0 时 的 系统 输出 幅度 比 4-0. 3 时 的 输出 
幅 值 下 降 约 20%。 可 见 ， 乘 性 方 波 信号 与 加 性 方 波 信号 对 信 噪 比 的 影响 是 不 同 的 。 实际 上 ， 
文 [13-14] 的 研究 结果 表明 ， 较 高 的 加 性 方 波 信和 号 幅 值 可 以 提升 系统 信 噪 比 , 而 由 图 2 和 表 2 
可 见 ， 较 高 的 乘 性 方 波 信和 号 幅 值 会 抑制 系统 输出 性 能 。 另 外 ， 信 噪 比 随 两 个 白 噪 声 间 的 耦合 
虽 度 A 变化 也 非 单调 变化 ， 如 图 2 和 表 2, “4A < 0.85, 信 噪 比 随 着 4 的 增 大 而 减 小 ， 而 
当 A > 0.9 时 ， 信 噪 比 随 着 4 的 增 大 而 增 大 。 


表 2 噪声 耦合 强度 4 和 波 信号 幅度 A 取 不 同 值 时 的 信 噪 比 (10 ) 
Table 2 The SNR for different values of coupling strength 4 between the noises and the amplitude A of the multiplicative 
square-wave signal 


A 0.30 | 0.50 | 0.65 | 0.70 | 0.75 | 0.80 | 0.85 | 0.90 0.95 1. 00 
A=0, 65 | 3.17 | 3.29 | 3.36 | 3.33 | 3.25 | 3.13 | 2.97 | 2.76 2.51 2. 26 
A=0. 67 | 2.95 | 3.07 | 3.19 | 3.20 | 3.17 | 3.10 | 2.98 | 2,81 2. 60 2. 36 
A=0. 69 | 2.76 | 2.86 | 3.02 | 3.05 | 3.06 | 3.03 | 2.96 | 2.83 2.67 2. 46 
A=0.72 | 2.51 | 2.58 | 2.77 | 2.83 | 2.88 | 2.90 | 2.88 | 2.83 2.72 2.08 
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图 2 噪声 耦合 强度 4 取 不 同 值 时 信 品 比 与 方 波 信号 幅度 A 的 关系 曲线 ， 参 数 为 a=1.9,b=0.7,D=0.35 , P=0.03,7r=0.4, 
=0.45 ,=1 
Fig.2 The SNR versus the amplitude A of the multiplicative square-wave signal for a=1.9,b=0.7,D=0.35,P=0.03,r=0.4, 


€=0.45,A, =1 for different values of the coupling strength A between the noises 


系统 参数 b 取 不 同 值 时 ， 信 品 比 与 直流 偏 值 7 的 关系 曲线 如 图 3， 由 图 3 可 见 每 条 曲线 
呈现 了 一 个 峰值 , 即 广义 的 随机 共振 现象 。 该 现象 与 文 [15] 中 的 结果 类 似 。 另 外 ， 容 易 看 
， 信 噪 比 也 是 参数 已 的 非 单调 函数 。 图 4 和 图 5 清晰 地 表明 ， 随 着 乘 性 噪声 强度 的 增 大 ， 
信 品 比 可 以 出 现 一 个 共振 峰 , 即 随机 共振 现象 , 该 现象 与 乘 性 噪声 作用 下 双 稳 系统 [12-14, 16] 
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中 出 现 的 现象 相同 。 而 且 ， 系 统 参数 4 非 单调 地 影响 信 品 比 , 如 图 4。 


图 3 系统 参数 b 取 不 同 值 时 信 噪 比 与 偏 值 > AMR, BRA a =2,D=0.8,P=0.04,24=04,€=0.8,A=0.01, A, =1 
Fig.3 The SNR versus the bias r for a=2,D=08,P=0.04,1=04,¢=0.8,A=0.01,4, =1 for different values of the 


system parameter b 


图 4 系统 参数 a 取 不 同 值 时 信 噪 比 与 乘 性 增 大 强度 D 的 关系 


线 , 参数 为 b=0.7 , P =0.25,A=0.15,r=0.1,s=0.1, A=0.1, 


为 =1 


Fig.4 The SNR versus the multiplicative white noise intensity D for p=0.7,P=0.25,4=0.15,r=0.1,s=0.1,4=01,4=1 


for different values of the system parameter a 


SNR 


图 5 噪声 耦合 强度 4 取 不 同 值 时 信 噪 比 与 乘 性 白 噪声 强度 D 的 关系 曲线 ， 参 数 为 a=0.9 ,b=0.9 ,P=0.2,r=0.01 ,c=0.1， 
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Fig.5 The SNR versus the multiplicative white noise intensity D for a=0.9,b=0.9,P=0.2,r=0.01,¢=0.1,A=0.1,4, =1 


for different values of the coupling strength A between the noises 


图 6 表明 ， 信 噪 比 是 加 性 白 噪 声 强度 的 非 单调 函数 。 随 着 加 性 噪声 强度 的 增 大 ， 系 统 输 
出 信 噪 比 先 单调 增 大 ,直到 其 最 大 值 ， 而 后 随 着 加 性 噪声 强度 的 进一步 增 大 而 单调 减 小 。 故 
相对 于 无 噪声 情形 ， 适 量 的 加 性 噪声 强度 可 以 提升 系统 输出 信 噪 比 。 


p 


图 6 系统 参数 b 取 不 同 值 时 信 噪 比 与 加 性 噪声 强度 P 的 关系 曲线 ， 参 数 为 a=2,D=0.2,4=0.3,r=0.2,g=0.2,A=0.1， 
Ay =2 


Fig.6 The SNR versus the additive white noise strength P for a=2,D=0.2,2=0.3,r=0.2,¢=0.2,A=0.1,4, =2 for 


different values of the system parameter b 


Sek, MERED Bea SWRA, (aR Cy DAB “SKE “a MEL (对 于 
A=0.72 情形 , 在 Æ0.8 时 出 现 最 大 值 为 2. 9， 生 0. 35 时 出 现 最 小 值 2. 5) ， 而 随 着 其 他 参 
数 的 变化 仅 出 现 一 个 极 值 ， 即 最 大 值 。 对 于 加 性 噪声 情形 ,在 信 噪 比 与 双 值 噪声 幅度 的 曲线 
上 或 者 仅 有 一 个 极 值 ， 即 最 大 值 ， 或 者 虽然 有 两 个 极 值 出 现 ， 但 先 出 现 最 小 值 ， 而 后 出 现 最 
大 值 。 另 外 ， 以 前 的 研究 表明 ， 较 高 的 加 性 方 波 信 号 幅 值 可 以 提升 系统 信 品 比 ， 而 本 文 的 研 
完结 果 发 现 ， 较 大 的 乘 性 方 波 信号 幅度 会 抑制 系统 输出 信号 ， 对 于 4=0.65 情形 ， 和 1.0 
时 系统 输出 幅度 比 A=0. 3 时 的 输出 幅 值 下 降 约 29%; 反之 ， 较 大 的 乘 性 双 值 噪声 幅度 则 可 以 
提升 系统 输出 信 品 比 ， 对 于 23.8 IB, €=3.5 时 系统 输出 幅度 为 E=0.25 时 的 输出 幅 
值 的 2.5 倍 。 


3 总 结 

本 文 研 究 了 乘 性 方 波 信 号 和 乘 性 双 值 噪声 作用 下 的 双 稳 系统 ， 观 察 到 了 两 种 随机 共振 现 
象 , 即 传统 的 随机 共振 与 广义 的 随机 共振 。 当 信 噪 比 随 着 噪声 强度 变化 时 ， 包 括 乘 性 双 值 噪 
声 强 度 、 乘 性 白 噪 声 强度 和 加 性 白 噪 声 强度 ， 乘 性 传统 的 随机 共振 现象 。 广 义 的 随机 共振 出 
现在 信 噪 比 与 系统 参数 、 直 流 偏 置 及 乘 性 方 波 信号 幅度 的 关系 曲线 上 。 那么 进一步 将 研究 结 
果 推 广 , 研究 拟 在 天 文 信号 处 理 方面 引入 本 研究 分 析 , 例如 在 不 同 接收 机 性 能 、 无 线 电 干扰 
膏 号 注入 情况 下 ， 对 接收 机 噪声 的 影响 情况 等 。 
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Abstract: Square-wave signal and dichotomous noise are widely used in scientific fields, and multiplicative signal is widely used in 
engineering practice. However, little study has been done on the nonlinear influence of multiplicative square-wave signal and 
multiplicative dichotomous noise on a bistable system, so in this paper we investigate the stochastic resonance phenomenon in a bistable 
system driven by multiplicative square-wave signal and multiplicative dichotomous noise is investigated. Under the adiabatic condition, 
the output signal-to-noise ratio (SNR) for the system is obtained. It is found that the SNR can obtain one maximum value and one 
minimum value with increasing the amplitude of the multiplicative square-wave signal and multiplicative dichotomous noise. Large 
values of the amplitude of the multiplicative square-wave signal can depress the system SNR, in the case of noise correlation rate 
A=0.65 , the system output amplitude for square-wave amplitude A=1.0 declines almost 29% by comparison with that for A=0.3, while 
a more amount of the multiplicative dichotomous will improve the system output, in the case of system parameter b=3.8, the system 
output amplitude for multiplicative dichotomous noise intensity ¢=3.5 is 2.5 times that for ¢=0.25 . One resonance peak can be 
observed when the system SNR varies with the multiplicative white noise, with the additive white noise, as well as with the system 
parameters. 


Key Words: Stochastic resonance; bistable system; multiplicative square-wave signal; multiplicative dichotomous noise 


